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The distribution of gluon fields in a baryon is of fundamental interest in QCD. We have observed 
the formation of gluon flux-tubes within baryons using lattice QCD techniques. In particular we 
use a high- statistics approach, based on the translational and rotational symmetries of the four- 
dimensional lattice, which enables us to observe correlations between the vacuum action density 
and the quark positions in a completely gauge independent manner. This contrasts with earlier 
studies which needed to use gauge-dependent smoothing techniques. We use 200 C(a 2 )-improved 
quenched QCD gauge-field configurations on a I6 3 x 32 lattice with a lattice spacing of 0.123 fm. 
Vacuum field fluctuations are observed to be suppressed in the presence of static quarks such that 
flux tubes represent the suppression of gluon-field fluctuations. We considered numerous different 
link paths in the creation of the static quark sources in order to investigate the dependence of the 
flux tubes on the source shape. We have analyzed 11 £ shapes and 8 Tand Y shapes of varying sizes 
in order to explore a variety of flux-tube topologies, including the ground state. At large separations, 
Y-shape flux-tube formation is observed. T-shaped paths are observed to relax towards a Y-shaped 
topology, whereas L-shaped paths give rise to a large potential energy. We do not find any evidence 
for the formation of a A-shaped flux-tube (empty triangle) distribution. However, at small quark 
separations, we do observe an expulsion of gluon-field fluctuations in the shape of a filled triangle 
with maximal expulsion at the centre of the triangle. Having identified the precise geometry of 
the flux distribution, we are able to perform a quantitative comparison between the length of the 
flux-tube and the associated static quark potential. For every source configuration considered we 
find a universal string tension, and conclude that, for large quark separations, the ground state 
potential is that which minimizes the length of the flux-tube. The characteristic flux tube radius 
of the baryonic ground state potential is found to be 0.38 ± 0.03 fm, with vacuum fluctuations 
suppressed by 7.2 ± 0.6%. The node connecting the flux tubes is 25% larger at 0.47 ± 0.02 fm with 
a larger suppression of the vacuum by 8.1 ± 0.7% 

PACS numbers: 12.38.Gc, 12. 38. Aw, 14.70.Dj 



I. INTRODUCTION 

Recently there has been renewed interest in studying 
the distribution of quark and gluon fields in the three- 
quark static-baryon system. While the earliest studies 
were inconclusive [l| , improved computing resources and 
analysis techniques now make it possible to study this 
system in a quantitative manner Hi, HQ. In partic- 
ular, it is possible to directly compute the gluon field 
distribution @, H, [§] using lattice QCD techniques sim- 
ilar to those pioneered in mesonic static-quark systems 

[Emm. 

Similar to Okiharu and Woloshyn [8| , our first interest 
is to test the static-quark source-shape dependence of the 
observed flux distribution, as represented by correlations 
between the quark positions and the action or topolog- 
ical charge density of the gauge fields. To this end, we 
choose three different ways of connecting the gauge link 
paths required to create gauge-invariant Wilson loops as 
in our preliminary study 9]. In the first case, quarks are 
connected along a T-shape path, while in the second case 
an L-shape is considered. Finally symmetric link paths 
approximating a Y-shape such that the quark positions 



approximate the vertices of an equilateral triangle are 
considered. The latter is particularly interesting as the 
probability of observing a A-shape flux-tube (empty tri- 
angle) distribution is maximised in this equidistant case. 
The most recent lattice QCD studies are in favor of the 
Y ansatz when it comes to the confinement potential |2| 
at long distances [!, 0, Q and flux tube formation [||. 
We also note that a theoretical calculation using center- 
vortex QCD with gauge group 5/7(3) also favors the Y 
ansatz [7| . In order to better compare T and Y shapes we 
study configurations in which the quarks are in the same 
location and only the connection point changes. Such a 
direct comparison is not possible with the L-shape. 

Because the signal decreases in an exponential fashion 
with the size of the loop, it is essential to use a method to 
enhance overlap with the ground state of interest. We use 
strict three-dimensional APE smearing [l3| to enhance 
the overlap of our spatial-link paths with the ground- 
state static-quark baryon potential. Time-oriented links 
remain untouched to preserve the correct static quark 
potential at all separations. We present results for two 
levels of source smearing including a minimal smearing of 
10 APE steps and an optimal smearing of 30 APE steps. 
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The comparison of two levels of smearing reveals inter- 
esting insight into the manner in which the flux tubes 
evolve to become the ground state. 

Complete details of the construction of Wilson loops 
for the determination of the static quark potential and as- 
sociated gluon-field distributions are presented in Sec.HTl 
Sec. IIIII presents the three-point correlation function cal- 
culated to reveal the response of the vacuum to the pres- 
ence of quarks. Of particular importance is the insensi- 
tivity to the choice of metric such that the observation 
of vacuum field suppression is unambiguous. 

Excellent signal to noise is achieved via a high- 
statistics approach based on translational symmetry of 
the four-dimensional lattice volume and rotational sym- 
metries of the lattice as described in detail in Sec. IIVI 
This approach contrasts with previous investigations, 
which used gauge fixing followed by projection to smooth 
the links and resolve a signal in the flux distribution 

Sec. [V] presents the simulation results for sources con- 
structed from 10-sweep smeared links. We find substan- 
tial correlation between the source and the observed field 
distribution. While these sources provide information 
vital to linking the flux-tube distribution to the static 
quark potential, they are not adequate to isolate the 
ground-state potential at large quark separations. There- 
fore in Sec. I VII we complement these results with sources 
constructed from optimal 30-sweep smeared links. Here 
clear deformation of the T-shape source to a Y-shape 
configuration is illustrated. A quantitative analysis com- 
paring the static quark potentials is used to identify the 
ground state and in Sec. IVIII the position, radius and 
expulsion extent of the flux tubes and the node connect- 
ing them are determined. A summary of our findings is 
provided in Sec. IVIIII 



II. WILSON LOOPS 

To study the flux distribution in baryons on the lat- 
tice, one begins with the standard approach of connecting 
static quark propagators by spatial-link paths in a gauge 
invariant manner. APE-smeared spatial-link paths prop- 
agate the quarks from a common origin to their spatial 
positions as illustrated in Fig. [TJ The static quark prop- 
agators are constructed from time directed link products 
at fixed spatial coordinate, YliUt(x,ti), using the un- 
touched "thin" links of the gauge configuration. 

The smearing procedure replaces a spatial link, U^{x), 
with a sum of 1 — a times the original link plus a/4 times 
its four spatially oriented staples, followed by projection 
back to SU(3). We select the special- unitary matrix 
which maximises KetriU^Uj}), where is the 
smeared link, by iterating over the three diagonal SU (2) 
subgroups of SU(3) repeatedly. In practice, eight iter- 
ations over the three subgroups optimizes the projected 
link to about 1 part in 10 s . 

We initially repeated the combined procedure of 
smearing and projection 10 times, with a — 0.7. How- 




FIG. 1: Gauge-link paths or "staples," Ui, Un and C/3, forming 
a three-quark Wilson loop with the quarks located at fi, r*2 



and r*3. e and e 



a'b' c 



denote colour anti-symmetrisation at 



the source and sink respectively, while r indicates evolution 
of the three-quark system in Euclidean time. 



ever, it will become clear that 10 times is not enough to 
isolate the ground state of the three-quark static-quark 
potential at large quark separations. We have found that 
repeating the smearing-projection procedure 30 times is 
effective in isolating the ground-state. In this case the 
string-tension of the three-quark potential matches that 
of the quark-antiquark potential. We will present for 
comparison, results for both 10 and 30 steps of APE- 
source smearing. 

Untouched links in the time direction propagate the 
spatially separated quarks through Euclidean time. In 
principle, the ground state is isolated after sufficient time 
evolution. Finally smeared-link spatial paths propagate 
the quarks back to the common spatial origin. 

The three-quark Wilson loop is defined as: 



tj/ ^abc^a'b'c' jjaa' jjbb t - 

Q = 3! 



3 > 



(1) 



where Uj is a staple made of path-ordered link variables 



(2) 



Uj = Pexp ig / dx li A' 1 (x) 



and Tj is the path along a given staple as shown in Fig.Q] 
In this study, we consider several two-dimensional 
spatial-link paths. Figures [2] and [3] give the projection 
of T and L shape paths in the x-y plane. Figure H shows 
the paths considered in constructing the Y-shape. In all 
cases the three quarks are created at the origin, O (white 
bubble), then are propagated to the positions Q±, Q2 or 
Q3 (black circle) before being propagated through time 
and finally back to a sink at the same spatial location as 
the source (O). 
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FIG. 2: Projection of the T-shape path on the x-y plane. 
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FIG. 3: Projection of the L-shape path on the x-y plane. 




TABLE I: (x,y) coordinates for the three quarks considered 
for Y- and T-shape paths. We also show the coordinates of 
the Fermat point, F, and the average distance of the quarks 
from this point (r a ). Also shown is the average inter-quark 
distance {d qq ). 





(x, y) Coordinates (lattice units) 


Distance (fm) 


# 


Qi 


Qi 


Qa 


F 






1 


(1,0) 


(-1,1) 


(-1,-1) 


(-0.42,0) 


0.15 


0.27 


2 


(2,0) 


(-1,2) 


(-1,-2) 


(0.15,0) 


0.27 


0.46 


3 


(3,0) 


(-1,2) 


(-1,-2) 


(0.15,0) 


0.31 


0.53 


4 


(3, 0) 


(-2,3) 


(-2,-3) 


(-0.27,0) 


0.42 


0.72 


5 


(4,0) 


(-3,4) 


(-3,-4) 


(-0.69,0) 


0.57 


0.99 


6 


(5,0) 


(-4,5) 


(-4,-5) 


(-1.11,0) 


0.72 


1.25 


7 


(7,0) 


(-4,6) 


(-4,-6) 


(-0.54,0) 


0.88 


1.52 


8 


(8,0) 


("4, 7) 


(-4, -7) 


(0.04,0) 


0.99 


1.71 



FIG. 4: Projection of the Y-shape path on the x-y plane. 




FIG. 5: The diagonal path (solid line) is taken as the average 
of the two possible paths going around the box enclosing the 
end points. 



For the Y-shape we create elementary diagonal "links" 
in the form of boxes as shown in Fig. \5\ The lxl and 
1x2 boxes are the average of the two path-ordered link 
variables going from one corner to the diagonally oppo- 
site one. Taking both of these paths better maintains 
the symmetry of the ground state potential and there- 
fore provides improved overlap with the ground state. 
We also consider 2x3 boxes which are the averages of 
the possible paths connecting two opposite corners using 
lxl and 1x2 boxes. We may create further link paths in 
the future for use in bigger loops if it becomes necessary. 
Hence a diagonal staple is, in fact, an average of several 
"squared path" staples connecting the same end points. 

The quark coordinates considered for Y and T shapes 
paths are summarised in Tabic |TJ The Fermat point - 
the point within an acute triangle that minimizes the 
sum of the distances to the vertices - is also indicated. 
We note that the origin (0, 0) is not at the centre of these 
coordinates. Rather the coordinates are selected to place 
the quarks at approximately equal distances from each 
other. As one can see, the Fermat point, F, (also called 
the Steiner point in this particular case) is very close to 
the origin in most cases. In the table we also indicate 
the average inter-quark distance (d qq ) and the average 
distance to the Fermat point (r s ). 

Note that in comparison to the measures considered 
in Ref. [1], we have Ly = 3{r s ) and La = 3{d qq ). If 



the quark positions were on the vertices of an equilateral 
triangle, Ly and La would be linked in a simple linear 
fashion. In the quark configurations of Table [H the quan- 
tity (r s ) is always within one part in a thousand of the 
radius of the circle circumscribing the triangle, therefore 
(r s ) is a good measure of the size of the baryon. 

For the L-shape we have studied a small sample of 
11 configurations, with 10-step smearing. The following 
coordinates for the quarks are considered: Q\ = (0,£), 
Q2 = (0,0) and Qz — (^,0) where I goes from 1 to 
11. Because all L-shape triangles are simply a scaled 
version of a unique triangle, all distances in this shape 
are a given constant times the scale factor (here "£") 
as in an equilateral triangle. In lattice units, we have 
F = (*=&l,*=&t), (O = ^> (d gg ) = ^ and 
the radius of the circumscribed circle is . The L-shapes 
are far further away from the ideal equilateral triangle 
than the T and Y shapes and (r s ) is not a good measure 
of the size of the baryon. 



III. GLUON FIELD CORRELATION 

In this investigation we characterise the gluon-ficld 
fluctuations by the gauge-invariant action density S(y, t) 
observed at spatial coordinate y and Euclidean time t 
measured relative to the origin of the three-quark Wil- 
son loop. We calculate the action density using the 
highly-improved 0(a 4 ) three-loop improved lattice field- 
strength tensor [14| on four-sweep APE-smeared gauge 
links. While the use of this highly-improved action will 
suppress correlations immediately next to the quark po- 
sitions, it will assist tremendously in revealing the flux- 
tube correlations which are the central focus of this in- 
vestigation. 

Defining the quark positions as f\ , f% and r% relative to 
the origin of the three-quark Wilson loop, and denoting 
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the Euclidean time extent of the loop by r, we evaluate 
the following correlation function 

C(y;fi,f 2) f 3 ;T) = 

(W 3Q (f u f 2 ,f 3 ;r)S(y,r/2)) 
(W 3Q (n,f 2 ,f 3 ;r))(S(y,T/2)) ' 

where (• • ■ ) denotes averaging over configurations and 
lattice symmetries as described below. This formula cor- 
relates the quark positions via the three-quark Wilson 
loop with the gauge-field action in a gauge invariant man- 
ner. For fixed quark positions and Euclidean time, C is a 
scalar field in three dimensions. For values of y well away 
from the quark positions ri, there are no correlations and 
C-*l. 

This measure has the advantage of being positive defi- 
nite, eliminating any sign ambiguity on whether vacuum 
field fluctuations are enhanced or suppressed in the pres- 
ence of static quarks. We find that C is generally less 
than 1, signaling the expulsion of vacuum fluctuations 
from the interior of heavy-quark hadrons. 



IV. STATISTICS 

In this investigation we consider 200 quenched QCD 
gauge-field configurations created with the 0(a 2 ^mean- 
field improved Luscher-Weisz plaquette plus rectangle 
gauge action [15( on 16 3 x 32 lattices at f3 = 4.60. The 
long dimension is taken as being the x direction mak- 
ing the spatial volume 16 2 x 32. Using a physical string 
tension of a = (0.440 GeV) 2 = 0.981 GeV/fm, the qq 
potential sets the lattice spacing to a = 0.123(2)fm. In 
our previous study @ we used 12 3 x 24 lattices with the 
same action and parameter j3. 

To improve the statistics of the simulation we use var- 
ious symmetries of the lattice. First, we make use of 
translational invariance by computing the correlation on 
every node of the lattice, averaging the results over the 
four- volume. 

To further improve the statistics, we use reflection sym- 
metries. Through reflection on the plane x = we can 
double the number of T and Y-shaped Wilson loops. By 
using reflections on both the plane x = and y = we 
can quadruple the number of L-shaped loops. 

We finally use 90° rotational symmetry about the x- 
axis to double the number of Wilson loops. This means 
we are using both the x-y and x-z planes as the planes 
containing the quarks. 

In summary, the Y and T shape Wilson loop results 
are the average of 16 3 x 32 x 2 x 2 x 200 = 104, 857, 600 
Wilson loop calculations. 



V. SIMULATION RESULTS WITH 10 SWEEPS 
SMEARING 

A. Qualitative results 

We begin by examining the correlation function 
C(y; fi, r 2 , rg; r) for 10-sweep smeared T-, Y- and L- 
shape sources at r = 2. Larger time evolutions display 
the same features and are discussed further below. 

At small quark separations, both the T- and Y-shape 
sources reveal a similar response of the vacuum to the 
quark positions. As seen in Figs. [5] and El isosurfaces 
of equal C(y) form a bubble for quark positions in the 
third configuration of Table HI Similar results are seen 
for smaller separations. 

The isosurface for lower levels of field expulsion (i.e. 
further away from the quark positions) displays an ap- 
proximately spherical shape as suggested by the surface 
plot for values of y just outside the quark positions. 

In these figures we also show the value of C (y) for y in 
the quark plane, (yi,j/2,0), on a surface plot (or rubber 
sheet). We clearly see that C(y) is less than 1 in the 
bubble, indicating an expulsion or suppression of vacuum 
gluon-field fluctuations. 

These features are also seen to some extent for the L 
shape up until £ = 3 as shown in Fig. [8] These results 
correspond to a radius of about 0.3 fm or under. 

At these small separations, one has the luxury of ex- 
ploring the Euclidean time evolution of C(y; ?\, r*2, r 3 ; r). 
While the depth of the vacuum field expulsion increases 
in going from r = 2 to 4, it saturates in going from 4 to 
6. In all cases however, there is no discernible change in 
the shape of the expulsion. Normalized contour plots of 
C(j7; t) are independent of r to an excellent approxima- 
tion. 

With the fourth set of quark coordinates the bubble 
loses more of its symmetry and a filled triangular shape 
emerges as in Fig. HI which corresponds to a radius of 
about 0.4 fm. This transition occurs at a radius of 0.35 
fm for the L shape. However, this shape is much less sym- 
metric to start with, as one of the inter-quark distances 
is 41% larger than the other two. 

It is interesting to note the expulsion indicated by the 
suppression of C(y) is largest at the centre of the quark 
configuration for all four quark-source shapes illustrated 
in Figs.[6]through[9] As such, there is no evidence of a A- 
shape flux-tube (empty triangle) distribution forming at 
small distances. Rather, the vacuum expulsion observed 
in Fig. [5] is shaped as a solid or filled triangle as opposed 
to tubes. Here the inter-quark distance is more than 0.7 
fm, which is outside the region where the potential is 
dominated by one-gluon exchange. 

For radii larger than 0.5 fm we clearly see the forma- 
tion of flux tubes as shown in Figs. [TU] to [T3] Here, a 
robust dependence of the observed flux-tube shape on 
the shape of the source is revealed. However, such a cor- 
relation is not surprising as these figures correspond to 
t = 2. In principle, with sufficient time evolution these 
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FIG. 6: Expulsion of gluon-field fluctuations from the region 
of static quark sources illustrated by the spheres. An iso- 
surface of C(y) is illustrated by the translucent surface. A 
surface plot (or rubber sheet) describes the values of C(y) for 
y in the quark plane, (j/i, 2/2, 0). Results are for a 10-sweep 
smeared T-shape source with quark positions as in the third 
configuration of Table[l] The maximum expulsion is 9.7% and 
the isosurface is set to 4.9%. 




FIG. 7: Isosurface and surface plot of C(y) for a 10-sweep 
smeared Y-shape source with quark positions as in the third 
configuration of Table[I] The maximum expulsion is 9.8% and 
the isosurface is set to 4.9%. Further details are described in 
the caption of Fig. [6] 




FIG. 8: Isosurface and surface plot of C(y) for a 10-sweep 
smeared L-shape source with quark separations of £ = 3. The 
maximum expulsion is 9.5% and the isosurface is set to 5.1%. 

ta,~j-i, — ^ ;u„,j ;~ — — ~c T7:„ Icl 




FIG. 9: Isosurface and surface plot of C(y) for a 10-sweep 
smeared Y-shape source with quark positions as in the fourth 
configuration of TableU] The maximum expulsion is 8.9% and 
the isosurface is set to 4.5%. Further details are described in 
the caption of Fig. [6] 



images will relax into the same ground-state configura- 
tion. In practice this is not possible as the signal is lost 
rapidly. Instead, one turns one's attention to improving 
the source to best isolate the ground state, as we shall 
do in the following sections. 

Nevertheless, some evolution of the flux tubes away 
from the source shape towards the Y-shape is already 
apparent. For example the flux-tube node of the T-shape 
source has moved by at least one lattice unit towards 
the Fermat point as illustrated in the contour plot of 
Fig. Q]|] Some curvature in the "top" of the T (left-most 
in the figure) is also apparent. In contrast, the observed 
Y shape of Fig. [14] is quite close to the source shape and 
the observed node of the flux tubes, which is very close 
to the origin, appears to have moved towards the Fermat 
point. 

For the L shape results of Fig. [To] the issue is more 
subtle. While the node appears to have moved in the di- 
rection of the Fermat point to the next grid intersection, 
this observation might also be due to the suppression of 
the correlation near the quark positions associated with 
our use of a highly improved action. 

While, Figs. IrH through [Pol show the action-based cor- 
relation function results, we have also examined correla- 
tions of the Wilson loops with the topological charge, the 
squared colour magnetic field and squared colour electric 
fields, all of which are gauge invariant. The results are 
qualitatively similar and in some cases indistinguishable. 



B. Quantitative results 

We now move to a quantitative study of these results 
exploiting the knowledge of the node position from our 
qualitative analysis. We begin by determining the effec- 
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FIG. 10: Isosurface and surface plot of C(y) for a 10-sweep 
smeared T-shape source with quark positions as in the seventh 
configuration of Table[l] The maximum expulsion is 8.3% and 
the isosurface is set to 4.4%. Further details are described in 
the caption of Fig. [6] 




FIG. 11: Isosurface and surface plot of C(y) for a 10-sweep 
smeared Y-shape source with quark positions as in the seventh 
configuration of TableU] The maximum expulsion is 8.3% and 
the isosurface is set to 4.4%. Further details are described in 
the caption of Fig. [6] 




FIG. 12: Isosurface and surface plot of C(y) for a 10-sweep 
smeared L-shape source with quark separations of £ = 10. 
The maximum expulsion is 8.8% and the isosurface is set to 
4.4%. Further details are described in the caption of Fig. [S] 



tive three-quark potential for the various quark positions, 
source shapes and Euclidean time evolutions. The vac- 
uum expectation value for W^q is 

oo 

(W 3 q(t)) = ]T C n exp(-a V n r), (4) 

71=0 

where V n is the potential energy of the n-th excited state 
and C n describes the overlap of the source with the n- 
th state. The effective potential is extracted from the 
Wilson loop via the standard ratio 

If the ground state is indeed dominant, plotting V as a 
function of r will show a plateau and any curvature can 
be associated with excited state contributions. Statistical 
uncertainties are estimated via the jackknife method [i"6| . 

Our results for the various quark positions and source 
shapes are shown in Fig. Q1)] All small shapes are stable 
against noise over a long period of time evolution and 
even some of the largest shapes show some stability be- 
fore being lost into the noise. 

Robust plateaus are revealed for the first four quark 
positions of Table [T] for the T and Y shape sources. This 
suggests the ground state has been isolated and indeed 
the four lowest effective potentials of the T- and Y-shape 
sources agree. This result was foreseen in the qualita- 
tive analysis where Figs. [6] and [7] for the T- and Y-shape 
sources respectively displayed the same correlations be- 
tween the action density and the quark positions. 

Conversely, the disagreement between Figs. [TU1 and[TT1 
indicates the ground state has not been isolated in one 
or possibly both cases. Indeed the nontrivial slopes of 
the seventh effective potentials of Fig. [16] for the Y- and 
T-shape sources confirm this. On the other hand, the 
curves are sufficiently flat to estimate an effective poten- 
tial at small values of r, and given knowledge of the node 
position from our qualitative analysis, one can make con- 
tact with models for the effective potential. 

The expected f dependence of the baryonic potential 
is SI 

j<k J 

where Cf = 4/3, a is the string tension of the qq poten- 
tial and L is a length linking the quarks. There are two 
models which predominate the discussion of L; namely 
the A and Y ansatze. 

In the A-ansatz, the potential is expressed by a sum 
of two body potentials [J]. In this case L = La/2 = 
3(d qq )/2 where La is the sum of the inter-quark dis- 
tances. In the Y-ansatz [1, Q, L = Ly = 3(r s ) is the 
sum of the distances of the quarks to the Fermat point. 
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FIG. 13: Contour plot of C(y) for y in the quark plane, 
(yiiJ/2, 0), for a 10-sweep smeared T-shape source with quark 
positions as described in the seventh configuration of Table U 
The white diamond shows the position of the Fermat point. 




FIG. 14: Contour plot of C(y) for y in the quark plane, 
(yi , 2/2, 0), for a 10-sweep smeared Y-shape source with quark 
positions as described in the seventh configuration of Table U 
The white diamond shows the position of the Fermat point. 
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FIG. 15: Contour plot of C(y) for y in the quark plane, 
(j/i, j/2, 0), for a 10-sweep smeared L-shape source with quark 
separations of I — 10. The white diamond shows the position 
of the Fermat point. 




\Y shape] 
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T shape\ 



012345678 
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\L shape] 



012345678 
X 

FIG. 16: Effective static quark potential for 10-sweep smeared 
Y-, T- and L-shape sources. Some points are offset horizon- 
tally to increase readability. From bottom up, the lines cor- 
respond to shapes 1 through 8 of Table U 
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In Figs . [T7l and [T8l we show the potentials from the dif- 
ferent shapes at r = 1 as functions of the inter-quark 
separation (d qq ) and the average distance to the Fer- 
mat point (r s ) respectively. In the case of the L-shape 
source, it is impossible to distinguish between the two 
ansatz from these figures as there is a linear relation be- 
tween (dqq) and (r s ). However, this is not the case for 
the T- and Y- shape sources. A linear or almost linear 
behavior for the long-distance potential is observed for 
Fig. [T7] considering (d qq ) and similarly in Fig. [T51 using 
(r s ). However, neither of these methods for expressing 
the inter-quark separation are capable of collapsing the 
potentials arising from T- Y- or L- shape sources to a 
single curve. 

In Eq. [6j L is meant to be the value of a length of string 
connecting the quarks (times a geometric constant for the 
A ansatz). Figs. [131 to [T5l indicate that the actual node 
of the flux tube formation is not necessarily in the initial 
location shown in Figs. to Q] nor, in the case of the T 
and L shapes, at the Fermat point. Using our knowledge 
of the node position in the observed state (which may not 
have relaxed to the ground state), we can determine the 
length of string required to connect the quark locations 
to the observed node position. 

In the case of the Y-shape source, both the observed 
node and the source node are very close to the Fermat 
point. The difference in the total distance from the quark 
to the node is at most a few percent compared to Ly. In 
the case of the T-shape source the observed node has 
moved by one lattice unit from the top of the T towards 
the Fermat point. In the smallest shape this means the 
node is extremely close to the Fermat point and that it 
is in the same location as the Y-shape node. In the case 
of the largest T-shape source considered in Table [1] the 
node may have moved more than one lattice unit. 

In Fig. [Tni we plot the potential of the various source 
shapes and sizes as a function of the total length of string 
connecting the quarks to the observed node. The length 
of string is computed as follows: 

• For the T shape, it is computed as the total distance 
from the observed node moved one lattice unit from 
the top of the T towards the Fermat point. 

• For the Y shape, it is taken to be equal to Ly. 

• For the L shape, it is taken to be 2t for the corre- 
sponding shape. 

Using this measure on the x-axis, all the data points 
fall close to a single curve as displayed in Fig. \T§\ The 
three curves for T-, Y- and L-shape sources with the 
node movement included are extremely close to each 
other. The outstanding curve labeled "T (original node)" 
emphasizes the importance of including node movement 
identified by the three-point correlation function. This 
observation favours the Y-ansatz as providing the correct 
understanding of the nontrivial, nonperturbative dynam- 
ics of QCD. 




0.4 0.8 1.2 1.6 



FIG. 17: Effective potential at r = 1 for various 10-sweep- 
smeared source shapes plotted as a function of the inter-quark 
separation (d qq ) — La/3. 
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FIG. 18: Effective potential at r = 1 for various 10-sweep- 
smeared source shapes plotted as a function of the average 
distance to the Fermat point (r s ) = Ly/3. 
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FIG. 19: Effective potential at r = 1 for various 10-sweep- 
smeared source shapes plotted as a function of the length of 
string connecting the quarks to the node location. 
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The asymptotic part of all the curves follows a straight 
line as suggested by Eq. \5\ and we can extract a string 
tension from this part. From the slope of the asymp- 
totic part of these curves the string tension is a w 
1.11 ± 0.01 GeV/fm. This is somewhat different from 
the results from Takahashi et al. [2] of 0.89 GeV/fm and 
also from the value we used to fix the lattice spacing 
(0.981 GeV/fm). On the other hand, Ref. @ established 
that the Y-law describes the asymptotic potential with 
the same string tension as the q-q potential. 

Since the Y shape source is already in a configuration 
providing access to the shortest string length, we must 
conclude that other aspects of the flux tube distribution 
have yet to evolve to become the ground state. Aspects 
including the flux-tube radius and the extent to which 
vacuum field fluctuations are suppressed inside the flux 
tube are examined in the following sections. 



VI. SIMULATION RESULTS WITH 30 SWEEPS 
SMEARING 

Given that simple Euclidean time evolution to the 
ground state is difficult for large quark separations, we 
turn our attention to improving the source by adjusting 
the amount of spatial-link smearing used in the source 
construction. 

For example, Takahashi et al. [2| used a large amount 
of spatial smearing in their studies compared to our first 
study of 10 smearing-projection sweeps. We have found 
that an increase in the amount of smearing leads to a 
reduction in the value of the potential for large quark 
separations while the smaller separations (r s < 0.5 fm) 
remain unchanged. This demonstrates that we have in- 
deed isolated the ground state for the smaller shapes after 
10 smearing sweeps, but not for the larger shapes. 

We therefore investigated systematically the effect of 
the smearing on the string tension. We find that the 
string tension decreases for increased smearing until it 
reaches a plateau at 30 steps of smearing. The resulting 
value for a — 0.97±0.01 GeV/fm is in agreement with the 
value used to fix the lattice spacing, a = 0.98 GeV/fm. 

The increase in smearing provides excellent plateaus in 
the effective potential of the Y-shape source at the ear- 
liest times. This can be seen in Fig. I2U1 which provides a 
side-by-side comparison of 30 and 10 smearing sweeps for 
the effective potential obtained from the first three values 
of the time parameter r. Improved plateau behavior and 
lower effective potentials are observed for all separations. 

Similar results are observed for the T-shape source as 
illustrated in Fig. [2lJ Figure l22l shows the complete time 
evolution of the potential of Y and T shapes for the 30 
sweep smearing results. While it is clear that the ground 
state has not yet been reached for the largest T-shape 
source separations, it will be fascinating to observe how 
the T-shape source flux-tube distribution has evolved to 
produce a much lower effective potential. 

The gluon-field distributions for the smaller source 
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FIG. 20: Comparison of the potential of the Y shapes for 10 
and 30 steps smearing. From bottom up, the lines correspond 
to shapes 1 through 8 of Table U 
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FIG. 21: Comparison of the potential of the T shapes for 10 
and 30 steps smearing. From bottom up, the lines correspond 
to shapes 1 through 8 of Table U 



shapes are largely unchanged by the shift from 10 to 30 
steps of smearing and we do not reproduce them here. 
However interesting results are revealed in the larger 
source shapes. The seventh configuration of TableUis de- 
picted in Figs . I!?3l and [Ml for the T- and Y-shape sources 
respectively. 

The observed flux distribution of the T-shape source 
in Fig. [53] displays significant evolution towards the ob- 
served Y-shape of Fig. [M] (and not towards a A shape). 
The node connecting the quarks in the T-shape source 
has now moved two lattice units in the direction of the 
Fermat point rather than only a single unit as in the 
10-sweep smearing case. While this is not so obvious in 
Fig. [23] it can be clearly observed in the contour plot 
of Fig. [21] and directly compared to Fig. [T2J It is ap- 
proaching the observed Y-shape source result of Fig. [^Hl 
In contrast, the Y-shape source flux-tube distribution is 
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FIG. 22: Study of the stability of the effective potential over 
Euclidean time for 30-sweep smeared Y and T shapes. From 
bottom up, the lines correspond to shapes 1 through 8 of 
Table U 



qualitatively unchanged, as can be seen in Figs. [24] and 
[26] compared to Figs. [TT1and[T4l 

The isosurface of Fig. [24] and contour plot of Fig. [26] 
provide our best determination of the distribution of 
gluon flux tubes within the ground state of the static- 
quark baryon. 

Motivated by our earlier discovery and knowledge of 
the precise location of the flux-tube node, we plot the 
effective potential as a function of the total length of 
the piece of string connecting the quarks to the node in 
Fig. [37] The length of the piece of string is now computed 
as follows: 

• For the Y shape, it is taken to be equal to Ly . 

• For the T shape, we first compute the total distance 
from a node moved one lattice unit from the base 
of the T. 

• We also consider a curve for the T shape where the 
node is moved one lattice unit from the top of the T 




FIG. 23: Expulsion of gluon-field fluctuations from the region 
of static quark sources illustrated by the spheres. An isosur- 
face of C(y) is illustrated by the translucent surface. A surface 
plot (or rubber sheet) describes the values of C(y) for y in the 
quark plane, (i/i,j/2,0). Results are for a 30-sweep smeared 
T-shape source with quark positions as in the seventh config- 
uration of Table [T] The maximum expulsion is 6.4% and the 
isosurface is set to 4.4%. 




FIG. 24: Isosurface and surface plot of C(y) for a 30-sweep 
smeared Y-shape source with quark positions as in the seventh 
configuration of Table[T] The maximum expulsion is 6.4% and 
the isosurface is set to 4.4%. Further details are described in 
the caption of Fig. 1231 



for small shapes (1 to 5 in Table [J) and two lattice 
units for the larger shapes. 

These results which converge very well to a single curve 
are contrasted by a curve for the T shape where the 
length of string is computed from the source-node po- 
sition located at the top of the T (original node). 

As can be seen in Fig. [57] the latter effective potential 
differs from the effective potential of the Y shapes. On 
the other hand, the two curves for the T-shape source 
with the node moved one or two lattice steps are very 
close together with the curve for the Y shapes. From 
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FIG. 25: Contour plot of C(j7) for y in the quark plane, 
(j/i, 3/2, 0), for a 30-sweep smeared T-shape source with quark 
positions as described in the seventh configuration of Table U 
The white diamond shows the position of the Fermat point. 
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FIG. 26: Contour plot of C(y) for y in the quark plane, 
(y i,y2, 0), for a 30-sweep smeared Y-shape source with quark 
positions as described in the seventh configuration of Table U 
The white diamond shows the position of the Fermat point. 
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FIG. 27: Effective potential at r = 1 for various 30-sweep- 
smeared source shapes plotted as a function of the length of 
string connecting the quarks to the node location. 



this figure alone it is quite hard to decide which of the 
two moved T shapes is the closest to the Y shape. The 
following section providing a quantitative analysis of var- 
ious cross sections of C(y) will resolve this issue. 



VII. FLUX-TUBE PROPERTIES 

Here the flux-tube and node properties as revealed 
by the three-point correlation function analysis in 
C{y;r\,r2,r^;T) are quantitatively examined by consid- 
ering several cross sections of C . In all cases y is con- 
strained to the plane of the quarks (yi,J/2j0) an d it is 
convenient to refer to these components as (x, y, 0) where 
x is the coordinate of the long direction. For the remain- 
der of this section the origin of the coordinate system is 
placed at the original position of the source node con- 
necting the quarks in the computations. 

We begin with a study of C(x, y, 0; r) with r = 2 for 
a fixed value of x. Larger values of t are examined be- 
low. In Fig. [28] we show the cross-section of the observed 
flux tube for the seventh T-shape source. Each curve 
corresponds to a fixed value of x. At x — 2 we are at 
the edge of the observed node position. At x = 12 we 
are past the quark in the long direction which is located 
at x — 11. Results for both 10- and 30-sweep smeared 
T-shape sources are displayed. 

Similarly, in Fig. [221 we show the cross-section of a flux 
tube observed for the seventh Y-shape source. Again, 
each curve corresponds to a fixed value of x. At x = 
we are close to the position of the node. At x — 9 we are 
past the quark in the long direction which is located at 
x = 7 in this case. 

While the 10-sweep-smearing computations display 
stronger vacuum-field action suppression, the 30-sweep- 
smearing results display a larger flux-tube diameter. The 
latter also display a Gaussian shape in the vacuum-field 
suppression. We can also see enhanced action suppres- 
sion near the location of the node and just before the po- 
sition of the quark. While the latter feature is suppressed 
to some extent by our consideration of a highly-improved 
action density, the greater action suppression at the node 
is physically interesting. The convergence of the curves 
in the intermediate regime reflects the formation of a well 
defined flux tube. 

For the T shape, the tube is reasonably well defined for 
x = 4 through 7 in both the 10- and 30-sweep-smeared 
sources. For the Y shape similar features appear for 
x = 2 to 5, again for both source smearings. We empha- 
size that it is the 30-sweep Y-shape source that provides 
the best overlap with the ground-state potential, as it is 
this source that provides the smallest and flattest effec- 
tive potential. Therefore, the 30-sweep Y-shape source 
provides the best determination of the actual shape of 
the flux tube within the ground-state baryon system. 

We can estimate the radius of the flux-tube by fitting 

2/2 

a function of the form 1 — Ae~ v /r to selected curves 
in Figs. [55] and [5U For the T shape we choose to fit 



TABLE II: Values of the fit parameters of the function 1 — 
A exp(— y 2 /r 2 ) to the curves x = 5 of Fig. 1281 and x = 3 of 
Fig. 1291 r is reported in lattice units (LU) and fm. The last 
column indicates the area under C(y) = 1 for the fitted curve 
in units of relative-depth times LU. 



Shape 


A 


r (LU) 


r (fm) 


At^/tt 


T (10 steps) 


0.070(1) 


2.3(1) 


0.28(1) 


0.28(1) 


Y (10 steps) 


0.073(1) 


2.2(1) 


0.27(1) 


0.28(1) 


T (30 steps) 


0.052(1) 


3.3(1) 


0.41(1) 


0.31(1) 


Y (30 steps) 


0.055(1) 


3.1(1) 


0.38(1) 


0.30(1) 



x = 5 and similarly we select x = 3 for the Y shape. 
All the non-linear fits hereafter have been done using the 
Levenberg-Marquardt method [TtJ • The results are sum- 
marised in Table [ill The fits are very good as illustrated 
by Fig. EDI 

Our earlier observations about the depth and width of 
the flux-tubes for 10 versus 30 source-smearing sweeps 
are confirmed in Table [III There is a significant corre- 
lation between the smeared-source size and the observed 
flux-tube radius. Whereas the source radius increases by 
70% in going from 10 to 30 sweeps at a smearing fraction 
of 0.7, the observed flux-tube radius increases by only 
40% for the Y-shape source. This indicates that some 
Euclidean time evolution of the flux-tube diameter has 
occurred in the 10-sweep smeared case, but not enough 
to relax to the true ground state captured by the 30- 
sweep smeared Y-shape source. Our best estimate of the 
ground-state flux-tube radius from the 30-sweep Y-shape 
source is r = 0.38(3) fm, where the uncertainty includes 
both statistical errors and the consideration of neighbor- 
ing flux-tube cross-sections. Similarly, A = 0.055(2) for 
the relative vacuum expulsion in the flux tube of the sev- 
enth quark separation. 

Finally in the last column we have the cross-sectional 
area created by the tube, computed from the fit. De- 
spite rather different radii and depths observed for the 
two different smearing levels, this quantity is relatively 
uniform. The breadth of the 30-sweep-smeared sources 
does give rise to a larger expulsion of vacuum field fluctu- 
ations in the flux-tube of the action density distribution, 
despite the fact that these sources give rise to a lower 
static quark potential. 

The observed flux distribution of the 30-sweep T-shape 
source in Fig. [23] displays significant evolution towards 
the observed Y-shape of Fig. [2H (and not towards a A 
shape). Further proof of this node displacement is pre- 
sented in Fig. 1311 Here the asymmetry about x — 
provides an estimate of the statistical errors. In both the 
10- and 30-sweep cases on the x = line, we can see the 
response to the quarks on both sides of the top of the T 
at x = ±6 and less suppression of the vacuum action in 
the middle. This indicates the node has moved in both 
cases. However, in the case of 10 sweeps of smearing, 
the node is encountered already at x — 1. On the other 
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FIG. 28: Value of C(x, y, 0) for the seventh T-shape source 
of Table [I] as a function of y (the short direction of the quark 
plane) for various fixed values of x (the long direction). Re- 
sults from sources constructed from 10 (top) and 30 (bottom) 
sweeps of APE spatially-smeared links are compared. In these 
plots the origin is placed at the position of the node in the 
source. 

hand with 30 sweeps of smearing one observes a plateau 
which shows that the combined effect of the quarks and 
the node makes C(y) constant. At x — 2 where we pass 
through the position of the node, the dip from the node 
is clear but it is almost the same depth as the plateau at 

35=1. 

Another perspective on the node properties is obtained 
by plotting the value of C(x, y, 0) along the line where the 
node and one quark stand on the long axis, y — 0. We 
show this in Fig. [35] for several 30-sweep T- and Y-shape 
sources. In these plots the x = position indicates the 
original position of the node in the source. 

In these figures we can typically observe two minima: 
one for the node on the left and one inside the quark 
position on the right. We choose to plot only the shapes 
for which these two features are resolved. For smaller 
shapes the node and quark peaks are not separable. 

For the T shape the node moves first one lattice unit 
(shape 5) then two (shapes 6 to 8). On the other hand, 
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FIG. 29: Value of C(x,y,0) for the seventh Y-shape source 
of Table U as a function of y (the short direction of the quark 
plane) for various fixed values of x (the long direction). Re- 
sults from sources constructed from 10 (top) and 30 (bottom) 
sweeps of APE spatially-smeared links are compared. In these 
plots the origin is placed at the position of the node in the 
source. 



for the Y-shape sources the node stays at its original 
location, indicating that the Y-shape source is already 
accessing the ground state. 

A fit using the left-hand side of the nodes with an ex- 
ponential curve of the form 1 — A e~( x ~ x °> l r can provide 
a more precise measure of the size of the node and of its 
actual location. The results have a slight dependency on 
the location of the cut in the data. In Tabic Hill we show 
the results of the fit for the seventh 30-sweep T- and Y- 
shape sources for two cuts. In cut I, only data from —16 
up to the minima are kept. For cut II, data up to one 
lattice unit on the right of the minima are kept. Both 
fits are in good agreement with our data as is illustrated 
in Fig. El 

Our best estimate of the ground-state node radius from 
the 30-sweep Y-shape source is r = 0.47(2) fm, where 
the uncertainty includes both the statistical error and 
the systematic uncertainty associated with the data cut. 
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FIG. 30: Example fit of a Gaussian to the flux-tube cross- 
section as per Table ITT1 The data points are for the seventh Y 
shape with 30 steps smearing as in the lower part of Fig. [29] 
but only for x — 3. 
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FIG. 31: Value of C(x,y,0) for the seventh T-shape source 
of Table U as a function of y (the short direction of the quark 
plane) for various fixed values of x (the long direction) near 
the top of the T-shape source. Results from sources con- 
structed from 10 (left) and 30 (right) sweeps of APE spatially- 
smeared links are compared. The origin is placed at the po- 
sition of the node in the source. 



TABLE III: Values of the node fit parameters in the function 
1 — A exp(— (x — xq) 2 /r 2 ) constrained to the left-hand side of 
the curves for the seventh 30-sweep T- and Y-shape sources in 
Fig.|32]for two data cuts. Cut I considers data from x — — 16 
to the minima and cut II from —16 to one lattice unit on the 
right of the minima, r is reported in both lattice units (LU) 
and fm. 



Shape 


A 


xo (LU) 


r (LU) 


r (fm) 


T (Cut I) 


0.061(1) 


+1.6(1) 


3.7(1) 


0.45(1) 


T (Cut II) 


0.062(1) 


+1.8(1) 


3.9(1) 


0.48(1) 


Y (Cut I) 


0.062(1) 


-0.5(1) 


3.7(1) 


0.45(1) 


Y (Cut II) 


0.063(1) 


-0.2(1) 


4.0(1) 


0.49(1) 
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FIG. 32: Value of C(x, y, 0) for a variety of T-shape (top) 
and Y-shape (bottom) sources of Table U as a function of x 
(the long direction of the quark plane) for y = such that 
the centre of the node is encountered. Results are for the 30- 
sweep smeared sources. In these plots the origin is placed at 
the position of the node in the source. 



This value is 24(9)% larger than the flux-tube radius of 
r = 0.38(3) fm. In addition, we observe a 15(3)% increase 
in the amplitude of vacuum-action expulsion, A, in the 
node of the flux-tube relative to the expulsion observed 
in the flux tube itself, for the seventh quark separation. 
The amplitude of the relative action expulsion in the node 
is reasonably uniform for the various quark separations 
at A = 0.062(2). Here the uncertainty reflects statisti- 
cal and systematic uncertainties associated with the data 
cuts and various quark separations. 

In these results r is linked to the radius of the node 
which approaches 0.5 fm, a little larger than the flux-tube 
radius found earlier for the 30-sweep smeared sources, xq 
is indicative of the precise location of the node and Ta- 
ble |TTT] provides xo — 1.7(1) LU for the T shape, consis- 
tent with the results of Fig. [57] and in accord with the 
contour plot of Fig. [231 

For the Y shape, Table Hill shows the node is very close 
to the origin defined at the node of the source. The 
movement of the observed node position to the left is 
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FIG. 33: Example fit of a Gaussian to the flux-tube node as 
per Table IIIII The data points are for the seventh Y shape 
with 30 sweeps of smearing as in the lower part of Fig. [32] 

consistent with movement towards the location of the 
Fermat point. In this case the Fermat point is located 
at x = —0.54 and Xq for Cut I is in agreement with this 
value. 

On a final note, we explore the Euclidean time evo- 
lution of our results for the 30-sweep smeared Y-shape 
source, which we have established to access the ground 
state potential in a quantitative manner. The early 
plateau behavior displayed in the top plot of Fig. [22] for 
the static quark potential suggests that there should be 
little to no Euclidean time dependence in our observa- 
tions. Indeed this is the case for the shape of the observed 
vacuum action suppression. As emphasized in Sec. [V] an 
examination of normalized contour plots reveals no dis- 
cernible t dependence in the flux-tube distribution. 

However, the depth of the relative vacuum-action sup- 
pression does indeed display a r dependence which is as- 
sociated with the use of a highly-improved lattice action 
density composed of 4-sweep APE-smeared links. Fig- 
ure [34] displays curves similar to that reported in Fig. [32] 
but this time for the second through fifth quark separa- 
tions listed in Table [T] Results for r = 2, 4, and 6 are 
illustrated for shapes 2 and 3 while the signal quality al- 
lows us to consider only r = 2 and 4 for shapes 4 and 5. 
Here the Y-shape source is considered. 

As the curves of Fig. [34] pass directly through the cen- 
tre of the node and the flux-tube to the right of the node, 
the amplitude of the relative vacuum-action suppression, 
A, may be simply read off the plot. While there is sub- 
stantial enhancement of A in going from r = 2 to 4, the 
value for A saturates at r = 6. 

Hence, the effect of further Euclidean time evolution is 
to enhance the amplitude, A, of the relative vacuum- 
action expulsion. The enhancement factor is uniform 
over the various separations at 1.3(1). Applying this cor- 
rection to our earlier results for the seventh quark sepa- 
ration where the node and flux tube are well separated, 
we find A = 0.081(7) and 0.072(6) for the node and flux- 
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FIG. 34: Value of C(x,0,0) for various Y-shape sources of 
Table U where the Euclidean time dependence of the vacuum- 
action suppression can be examined. Here, x is along the 
long direction of the quark plane (z = 0) and with y — 
the centre of the node is encountered. Results are for the 
30-sweep smeared sources. 

tube respectively. 

VIII. CONCLUSIONS 

We have studied extensively correlations of the 
vacuum-action density with Wilson loops describing the 
three-quark static-baryon system as defined in Eq. ([3]). 
A high-statistics approach based on the translational and 
rotational symmetry of the four-dimensional lattice vol- 
ume is adopted to avoid the need for gauge-dependent 
smoothing techniques. 

The correlations indicate that the vacuum-action den- 
sity is suppressed in the interior and surrounding volume 
in accord with previous findings in the quark-antiquark 
case. "Flux tubes" represent the expulsion or suppres- 
sion of vacuum gluon-field fluctuations. 

For small quark separations, the ground state is ac- 
cessed easily through the use of spatially-smeared APE 
links in the construction of the source and sink of the 
three-quark system. As illustrated in Fig. [71 the vacuum- 
action suppression is quite spherical just outside the 
quark positions, and is strongest in the very centre of 
the three-quark system. As such, there is no evidence 
of the formation of a A-shape flux-tube (empty triangle) 
distribution. 

As the quarks become more separated the vacuum 
expulsion takes the shape of a solid or filled triangle. 
We conclude that at small quark separations i.e. for 
(r s ) < 0.5 fm Y-shape flux-tube formation does not oc- 
cur. As illustrated in Figs. 0E] and [9] a localized junction 



of flux tubes is not observed. 

For quark separations where the quarks are more than 
0.5 fm from the system centre, flux-tube formation is 
observed, in support of Refs. 0,0, Here, the action 

distribution shows some dependence on the source due 
to the limited Euclidean time evolution that is possible. 
Through a study of 8 T and Y shapes and 1 1 L shapes of 
varying sizes, we are able to quantitatively demonstrate 
that the observed potential is a linear function of the 
length of string required to connect the quark positions 
to the observed position of the node connecting the flux 
tubes. This is illustrated in Figs. fT9l and l27l for the 10- 
and 30-sweep sources respectively. As such, the ground 
state potential at large quark separations is the flux-tube 
distribution which minimizes this length. In this case the 
node will be located at the Fermat point. 

We have shown that T-shape sources have a flux-tube 
distribution which relaxes towards a Y-shape distribu- 
tion as the ground state of the static quark potential is 
accessed. This is illustrated in Fig. [55] with the corre- 
sponding effective potentials presented in Fig. [201 While 
the formation of a Y-shape flux-tube has been observed 
before (see @ for example), to the best of our knowledge 
this is the first time such a structure has been observed in 
a gauge invariant manner and observed to emerge from 
sources that do not have a Y-shape bias. 

We have successfully exposed the flux-tube distribu- 
tion of the ground-state static-quark baryon system at 
large separations. Figures [Ml and |2"61 provide the best il- 
lustrations of this flux-tube distribution. Animations for 
a variety of quark separations are available on the web 

m- 

Finally the characteristic sizes of the flux-tube and 
node have been quantified. We find the ground state 
flux-tube radius to be 0.38(3) fm with vacuum-field fluc- 
tuations suppressed by 7.2(6)%. The node connecting the 
flux tubes is 24(9)% larger at 0.47(2) fm with a 15(3)% 
larger suppression of the vacuum action at 8.1(7)%. 
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